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I study the average deformation rate of an amorphous material submitted to an external uniform 
shear strain rate, in the geometry known as the split-bottom configuration. The material is described 
using a stochastic model of plasticity at a mesoscopic scale. A shear band is observed to start at 
the split point at the bottom, and widen progressively towards the surface. In a two-dimensional 
geometry the average statistical properties of the shear band look similar to those of the directed 
polymer model. In particular, the surface width of the shear band is found to scale with the system 
height H as H a with a = 0.68 ± 0.02. In more realistic three dimensional simulations the exponent 
changes to a = 0.60 ± 0.02 and the bulk profile of the width of the shear band is closer to a quarter 
of circle, as it was observed to be the case in recent simulations of granular materials. 

PACS numbers: 



I. INTRODUCTION 

Amorphous and granular materials have a strong ten- 
dency to develop localized deformation when submitted 
to shear deformation. This deformation concentrates 
usually in narrow regions, called shear bands. Upon a 
uniform rate of shear deformation, the material expe- 
riences a finite, stationary deformation rate within the 
shear band, whereas the deformation rate is vanishingly 
small outside iti^ In some cases the existence of strain 
localization can be associated to a decreasing form of the 
stress as a function of strain rate in the material. This 
decreasing leads to instabilities that drive the formation 
of shear bands. 

An interesting limiting case corresponds to the so 
called quasi-static shear banding. In this case the ap- 
plied strain rate is vanishingly small, and the velocity 
field across the sample must change by a factor when 
there is a change of the driving force by the same factor. 
It seems that in this case, an stationary, finite width shear 
band can only be obtained through the application of 
appropriate conditions at the boundaries that break the 
uniformity of the system. In fact, simulations of homo- 
geneous amorphous materials well below the glass tem- 
perature show the existence of shear band which however 
thicken in time, eventually taking over the whole system^ 

A particular example in which the shear band is in- 
duced through the use of appropriate boundary condi- 
tions is the split-bottom container geometry introduced 
ini£ in the form of a modified Couette cell. In this con- 
figuration the bottom of the cell is split into two disks 
that rotate relatively to each other, and are attached to 
the inner and outer cylindrical wall of the cell. The sharp 
velocity interface at the bottom of the container induces 
an abrupt change of the velocity profile of the material 
close to this split line. Experiments using granular mate- 
rials show that in these conditions the material inside the 
cell reaches (on average) a stable configuration in which 
the transition between particles moving rigidly with one 
or the other part of the bottom surface occurs in a tran- 
sition zone that defines the shear band. The form of this 
zone is well defined and non-singular even in the limit 



of 'quasi-static' shear bands, i.e, when the relative veloc- 
ity of both parts of the container is vanishingly small. 
However, when horizontal planes progressively away of 
the bottom are considered (and thus the influence of the 
split-bottom decreases), the width of the shear zone be- 
comes larger, and no upper limit seems to exist. The 
width of the shear band starts as zero right at the bottom 
of the container, and progressively increases when mov- 
ing upwards, being maximum at the free surface. The 
detailed study of the form of this profile is the subject of 
this paper. 

The geometry that I will use here is a more symmetric 
version of the split-bottom experiment, which is depicted 
in Fig. [TJ The bottom of the container defines the y-z 
plane, which is split along the y — line. The x direction 
is perpendicular to the bottom plane. The system is in- 
finite in the z direction (in the three dimensional numer- 
ical implementation, periodic boundary conditions along 
z will be used). The configuration is symmetric upon a 
change of sign of the y coordinate, and thus the center of 
the shear band is always located at y = 0. Note that this 
applies to the time averaged position of the shear band, 
and not necessarily to an instantaneous configuration. 
We will come back to this fact in the discussions. This 
geometry has been analyzed numerically using frictional, 
spherical particles in£. 

To predict, or justify, the form of the shear band we 
have to make assumptions about the structural response 
of the medium. One of the simplest starting points to 
model a plastic material is the assumption of 'perfect 
plasticity'. Under this hypothesis, each piece of material 
behaves elastically at small deformation, i.e., the stress 
is proportional to the strain, but if some threshold de- 
formation is over passed, the stress remains constant. It 
can be shown that if this simple prescription is followed, 
a singular shear band is obtained that goes from the split 
line at the bottom to the free surface at the top J- If the 
threshold deformation is uniform in the sample the shear 
band is a vertical plane above the split line, namely the 
y-z plane. However, if different thresholds are assigned to 
different points of the sample the shear band obtained is 
the homogeneous projection in the z direction of a wavy 
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line in the x-y that goes between the split line at the bot- 
tom, and the free top surface. This shear surface can be 
obtained as the one that minimizes the energy generated 
as 'heat of friction' during sliding. 




FIG. 1: (a) The 'split bottom' geometry. Left and right parts 
of the system are driven at a very small velocity v and —v, 
respectively. The top surface is free of force. The velocity 
of the material match these values close to the bottom sur- 
face, and interpolates continuously in the bulk, as indicated 
schematically, (b) A 'split top-bottom' geometry. Now the 
top surface is also split into two, like the bottom. Note the 
definition of the height H in both cases, which is made in 
this way to allow a better comparison. Even in the absence 
of gravity effects is not obvious that the configuration in (a) 
should be equivalent to half of the configuration in (b), as it 
is explained in the text. 

In Ref.— , Torok et al. put forward the idea that this 
minimal path is not persistent in time, since upon some 
rearrangement of the particles, the thresholds near the 
shearing surface are likely to change, and the path can 
move around slightly. In this way, a fluctuating path 
connecting the split point and free surface is obtained. 
It was thus suggested that the average configuration of 
this fluctuating path produces the continuum shear band 
observed in the experiments. Simulations based on the 
two-dimensional version of this idea (in which the z co- 
ordinate is eliminated, and the problem is solved in the 
x-y plane, assuming homogeneity along z) find a shear 
band that, after time averaging, is an object that widens 
from bottom to top. 

An interesting magnitude that has been studied with 



some detail in this context is the scaling of the width 
w of the shear band at the surface as a function of the 
cell height H. In the first reported experimental worki it 
was found that there is a power law dependence for this 
quantity, with the shear band width w scaling as H a , 
with a between 0.5 and 1.0. In the numerical work of 
Torok et a value close to 2/3 was found for a , and 
it was noticed that this is the analytical value known to 
occur in the case of directed polymers (DPs) in 1 + 1 
dimensions 3 After this, it seems that people have taken 
the value of 2/3 as a reference, and systematically try to 
scale the results of experiments or simulations with this 
value. 

In my view there is not any strong reason why this 
should be the exact scaling exponent. The only point 
to be remarked is that the problem of DPs at zero 
temperature 9 is in essence very similar to the case of 
the fluctuating shear line in two dimensions. The main 
difference is that in the case of polymers the whole dis- 
tribution of thresholds (and not only close to the shear 
line) is updated in every time step. An additional differ- 
ence is that the model in£ is not 'directed'. However it 
is known that for polymers this restriction is not crucial 
either. The fluctuation of the tip of DPs at the surface of 
the system is analytically known to grow as H 2 /*. How- 
ever, the problem of the shear band is a continuous three 
dimensional one, and this makes the comparison with the 
DP rather obscure. 

Much less experimental information exists on the 
width of the shear band inside the material, namely w(x), 
for < x < H. This is in part originated in the diffi- 
culty to observe the movement of the particles in the 
bulk, although some advances have been made in this 
direction*^ In Ref£, Ries et al. did numerical simu- 
lations using a rather realistic three dimensional model 
of spherical particles interacting through friction contact 
forces. I consider that the results in this work comple- 
ment nicely the available experimental measurements. In 
these simulations, in addition to the (approximate) 2/3 
power scaling of w(H), a quarter of a circle profile was 
adjusted for the width in the interior of the material, 
namely w(x) ~ y / x(2H — x). 

It would be interesting to compare these results with 
the equivalent ones obtained with the previously men- 
tioned model of fluctuating shear lines^. Unfortunately, 
the results in Ref.— are not very detailed concerning the 
width across the whole sample. Coming to the problem 
of DPs, to my surprise I have not found a detailed study 
of the full profile either, and then I set up a simulation 
code to study this case too. The results obtained for DPs 
will be presented below. 

The DP model, or the model of fluctuating shear line 
of Ref. 8 consider that instantaneously, the shear band 
connects the split-bottom point and the free surface by a 
single line that fluctuates in time. Experimentally, there 
is no evidence that this is the case. In fact, even for 
the lowest velocities that can be measured, the velocity 
profiles that are observed seem to be continuous. This is 
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one reason to try some model which is a bit more realistic 
and considers the possibility of continuous deformations, 
both in time and space. Another reason to search for 
a more realistic model is the fact that the DP model, 
and the model of fluctuating shear line of Ref4 are emi- 
nently two dimensional, and cannot be easily generalized 
to three dimensions. 11 In fact, in three dimensions, an 
instantaneous slip surface cannot occur throughout the 
system, unless its profile is independent of the third (z) 
dimension. The problem of the width of a shear band 
we are discussing is really a three dimensional one. Al- 
though it can be investigated by two dimensional models 
(defined in the x-y plane in Fig. [I}, it is not clear if 
there are crucial issues associated to the third direction. 
In this sense, note that the problem is homogeneous in 
the z direction only on average, but not instantaneously 
due to the randomness of the material. The model pre- 
sented in the next Section can be implemented in three 
dimensions, and I will present results both for two and 
three dimensions. 



band of zero width. In order to obtain a shear band of 
finite width, the totally averaged, deterministic evolu- 
tion does not seem to be enough, a stochastic component 
must remain^ There are different possibilities for the 
choice of the stochastic evolution, but I think the kind 
of results that are obtained is rather independent of the 
details. The model adopted here is the following. I start 
with a zero local plastic strain. Upon evolution according 
to Eq. [2 the strain field changes, and also does the 
local energy / (Eq. [T]) . The plastic strain is fixed to zero 
until / overpasses some local threshold value fth(x,y), 
then the value of e pl (x,y) is set instantaneously to the 
local value of the strain. In this way the local free energy 
is relaxed to zero. In addition, a new, random value of 
fth{x,y) is chosen. This accommodation of the plastic 
strain can be considered to be the plastic events that oc- 
cur in the system, which allow to produce a potentially 
very large deformation, at the cost of a moderate energy 
consumption. 



II. A MESOSCOPIC SCALE MODEL 

The model studied here is a simplification over the 
model introduced in Ref.— to study the finite width of 
shear bands in the presence of structural relaxation. Here 
this structural relaxation will be absent. I consider a two 
dimensional planar geometry in 'mode III' configuration, 
in which an out-of-plane displacement u(x, y) is defined. 
The strain field is defined by the two-component vector 
ei = du/dx, e2 = du/dy. 

A perfectly elastic two dimensional material is modeled 
through the local free energy / = B(e^ + e^), where B is 
related to an elastic constant of the material. A plastic 
material possesses a plastic strain field e P \ (i = x, y), in 
such a way that the local free energy reads 

f(x,y) = B[(e x -eff + (e y -e^] 

F = J f(x,y)dxdy (1) 

The time evolution of u(x, y) can be considered to be 
determined in the quasi-static limit by a relaxational first 
order equation of the form 

du(x,y) = A SF 
dt 5u(x, y) 

The model must be complemented by the time evolu- 
tion of the plastic strain. In principle, a totally determin- 
istic model for the evolution of ef could be tried. For 
instance, we can postulate that if the local energy exceeds 
some threshold value, the system adjusts its plastic strain 
to reduce the local free energy. However, all attempts 
to implement this kind of deterministic model led in the 
present geometry to a singular slip surface, i.e., to a shear 



III. RESULTS 

A. Directed Polymers 

First of all I will present the results obtained using the 
DP model. As a first check, I reproduced with high preci- 
sion the analytical 2/3 power for the scaling of the width 
at the surface if the number of elements of the polymer is 
sufficiently large, as can be seen in Fig. [5J The full curve 
uu(x) of the largest polymer studied is also plotted in that 
figure. It coincides at the surface with the previous curve, 
by definition, but is systematically different in the bulk. 
To my knowledge, there is no analytical expression for 
the w(x) profile of a DP. Note that the half-profile for a 
DP fixed at both ends (plotted also in Fig. [5J) shows clear 
differences with the free end case: the width of the free 
end case is systematically larger than the fixed end case, 
they only coincide asymptotically close to the fixed point. 
Note also that the curve for the two fixed ends case has 
clearly zero slope at the middle point, which is obviously 
due to symmetry reasons, but this is not obviously the 
case for the free end case. 

To analyze these results in more detail, it turned out 
to be very successful to apply the theory of Gaussian sig- 
nals developed in Refi^. This consists in the study of 
the properties of a self-affine stochastic process, submit- 
ted to the appropriate boundary conditions, assuming a 
decomposition into uncorrelated Fourier modes of Gaus- 
sian amplitude. Although the hypothesis of the modes 
being uncorrelated is usually not satisfied, this analysis 
is known in particular cases to provide results that fit 
the exact data with high accuracy.— The simplest case 
to start with is that of the fixed ended polymer. In this 
case it is assumed that the stochastic process can be de- 
composed in Fourier modes in the form: 
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□ Surface width, free end (H=x) 
O Bulk width, free end 
A Bulk width, fixed end 



X 



w(x) 



V2H 



Vx(H ~ x) 



(6) 



i.e., the well known half circle, which is the exact result 
in this case. For DPs we have a = 2/3, and the result for 
w(x), fitting the global factor Co is shown on top of the 
numerical results in Fig. [5] and cannot be distinguished 
from the numerical data. 

Going to the case of the free end polymer, note that 
in the problem of thermally fluctuating elastic line, only 
modes that reach the surface with zero derivative must 
be considered. This leads to a profile of the form 



FIG. 2: (Color online) Squares: Numerical results for the 
average fluctuation of the free end of a directed polymer as a 
function of the number of its elemental pieces. The analytical 
2/3 power (thin continuous) is very well reproduced once the 
number of elemental pieces becomes large. The full profile of 
a free end polymer (circles) , and the half profile of a fixed end 
polymer (triangles) with H — 600 are also plotted (for better 
visualization not all points are plotted). Continuous lines: 
fitting using the random Gaussian signals of Ref.J^ (Eq. [9]). 



Y(x) = g n A n sm(mrx/H) 



(3) 



where Y(x) is the y coordinate of the polymer at height 
x (see Fig. 1), g n are uncorrelated random numbers with 
zero mean and unitary variance, H is the distance be- 
tween the two fixed ends of the polymer, and A n is the 
amplitude of the corresponding mode. For self affine pro- 
cesses with scaling properties, the amplitudes A n must 
scale as a power of n, namely A n — Con -7 , where an 
-//-dependent global constant was introduced. By multi- 
plying and taking averages over realizations we get the 
statistical width w at distance x of one fixed end as 



w 2 {x) = C, 



\n—l 



" 27 sin 2 



(rnrx/H) 



(4) 



where I have made use of the fact that different g n are 
uncorrelated. The exponent 7 can be related to our expo- 
nent a for the expected scaling of the stochastic process 
at the surface. It is not difficult to show that a = 7 — 1/2. 
Thus we can write 



\x) = C% [J2 n ~ 2Q_1 sin2 (n-Kx/H) 



(•5) 



\n=l 



As an example, note that the thermally fluctuating 
elastic line problem, has a = 1/2, the sum in Eq. (JSJ) 

gives 



w\x) = Cl 



E 



n-^sin 2 



(nirx/2H) (7) 



\n=l,3,5,... 



where now the sum is over half the modes we had before. 
The result obtained scales as 



w(x) 



(8) 



In particular, this profile does not have zero derivative at 
the surface, although it is constructed with modes that 
have zero derivative individually. 

Coming back now to the free end DP polymer, an at- 
tempt can be made to fit the curve of an open end poly- 
mer weighting appropriately the modes with zero value 
at the surface, and those with zero derivative, keeping 
the self-affine nature of both of them individually. This 
corresponds to a width function of the form: 



™ 2 (z) = c 2 (X,(i + /?(-i) n K 



sm 2 (mrx/H) 



The best fitting to the numerical results is obtained 
when we choose (3 = 1/2, and Co as the same value used 
to fit the fixed end case. The fitting, shown in Fig. [21 
cannot be distinguished from the numerical results either. 
Note that having the same value of Co means that the 
asymptotic form of w(x) for x — > is exactly the same 
for the free end and open end situations (this is because 
the alternating in sign (3 term in © introduces for small 
x a contribution that is a larger power of x than the main 
term). 

A very sensitive quantity that can be used to evaluate 
the results just presented is the local slope s of the w(x) 
curve, defined as s(x) = d\nw(x)/dhix. These results, 
both from numerics and using the fitting with Gaussian 
signals are shown in Fig. [3l Two facts must be no- 
ticed: first, for open end polymers, the slope vanishes at 
the surface as (H — x) 1 / 3 , more abruptly than the linear 
vanishing of the fixed end case at the mid point. Sec- 
ond, although the convergence is slow with system size, 
data seem to indicate that the form of w(x) for x — > 
is w(x) ~ x 2 / 3 , i.e., it has the same exponent that the 
surface scaling (which is, by the way, a property of the 



5 




x/H 



FIG. 3: (Color online) Logarithmic derivative of the bulk 
width data of free end DPs with H = 600 and 1200, and 
fixed end DP with H — 600, and the corresponding fitting 
according to random Gaussian signals (— ). Note that for free 
end DPs, the derivative vanishes at the surface at (H — - x) 1 ^ 3 . 
Note also that the data for small arguments tend to approach 
the trend of the continuous line as size is increased, although 
very slowly (the blow-up shows only data for the open end 
case) . 



Gaussian signal too). I will come back to the meaning of 
this fact in the last section. 

At present, I do not have any explanation that jus- 
tifies the accuracy of the fitting proposed, and neither 
an explanation for the remarkably value of (3 obtained 
in the fitting of the open end polymer. The results just 
presented for DPs will serve as a reference for comparing 
those obtained with the mesoscopic plastic model. 



B. Plastic mesoscopic model: two dimensions 

In a first stage, I implemented the plastic mesoscopic 
model on a two dimensional rectangular geometry corre- 
sponding to the x-y plane in Fig. [TJ in which the border 
of the right half is given some velocity v, and the border 
of the left half is given the opposite velocity —v. The 
value of v is set small enough for the total velocity field 
to be simply proportional to this value. 

For the split bottom geometry, the top surface is left 
free, i.e., a condition of zero normal stress is applied. In 
some results presented, the top surface will be also split 
in the same manner that the bottom surface, and this 
will be a 'split top-bottom' configuration. Note that al- 
though I use the expressions 'top' and 'bottom' to agree 
with the experimental device, in the present case I do not 
include a gravitational field that singles out a vertical di- 
rection. This choice is made as the numerical simulations 
in- indicate that gravity plays only a minor role in the 
structure of the shear band. Simulations were performed 
for different values of the total height of the system H, 
and in each case the size along the y direction was chosen 
sufficiently large for the shear band not to be affected by 



the lateral walls (typically, the size along y is around ten 
times the expected value of w at the surface) . 

Starting from rest, with both and e? 1 set to zero, the 
simulation is initiated by applying a finite velocity ±v to 
the left and right borders. At the first stages of the sim- 
ulation, the system is seen to respond elastically, since in 
no place the value of fth is over passed. At some moment, 
fth is over passed for the first time at some place, the cor- 
responding value of e? changes, and this feeds back to 
modify the evolution of e^. After many plastic events 
at different spatial positions, the system reaches a sta- 
tionary state, (on average, since stochastic fluctuations 
of the plastic strain remain), and at this stage statistics 
can be accumulated to get the average velocity field. In 
the simulations presented, the threshold values fth are 
chosen randomly in the interval 0.5/ t ° h < fth < l<5/$, 
where f® h is a reference value. Dimensionless units are 
set by taking f° h = 1 and XB = 1 in Equations HJ)-©- 
In these units, the small velocity used to drive the system 
is 10~ 4 , and the time step y 0.1. 

The shear band forms as an effect of the sequence of 
plastic events in which the values of ef Z are adapted at 
different points in space. These events are correlated in 
space, in general. A plastic slip occurring at some po- 
sition (x, y) reinforces the possibility of a similar event 
occurring above or below, at the same y coordinate, and 
displaced in the x direction. However, the stochastic na- 
ture of the model makes this overall tendency not to be 
strict. Accumulated plastic deformations over progres- 
sively larger periods of time, shown in Fig. 01 allow to 
understand the mechanics of the process. On very short 
time scales (Fig. BJa)), individual events at isolated po- 
sitions are observed. At intermediate time scales (Fig. 
HJb) and (c)), we see how these events are correlated 
in space to generate a sort of finite length shear lines. 
Note however, that contrary to the fluctuating shear line 
model, here the plastic deformation does not accumulate 
at a single shear line, irrespective of the time interval 
considered. For the largest time scales (Fig. Hljd)) we 
can see the tendency for the the fluctuating behavior to 
be averaged out, and for sufficiently large time intervals 
we get eventually a well defined, continuous shear band. 

Now I will focus on the properties of the averaged shear 
band in this stationary situation. In Fig. [5fa) I plot the 
full profile for the width w of the shear band obtained 
in systems of different heights H . These data can be 
very well scaled if plotted as w/H a against x/H. This is 
basically the same result obtained in experiments and in 
other simulations, extended here to the full profile inside 
the sample. The best fitting for the power a gives a = 
0.68 ± 0.02 (see Fig. Mb)). Applying to this case the 
fitting using Gaussian signals, we get the fitting shown 
in Fig. [SJ This is obtained using expression ^ with a = 
0.68 (taken from the surface scaling), and (3 = 0. The 
corresponding logarithmic derivative of the data and the 
fitting are shown also in Fig. [f^b). Although a slightly 
different from zero value of f3 gives a somewhat better 
fitting, I think it is within the numerical uncertainty. The 
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FIG. 4: Schematics of plastic deformation in the horizontal 
direction (du/dx) accumulated over progressively larger peri- 
ods of time At, namely At /to = 10 2 , 10 3 , 10 4 , and 1.5 x 10 5 , 
from (a) to (d). Darker regions indicate zones with larger ac- 
cumulated plastid deformation. Note that the plastic events 
appear at rather isolated positions (a), that then mersge to 
generate pieces of shear lines (b-c) . Over very large periods of 
time the fluctuating behaviors averages out, and a stationary, 
continuous shear band is observed (d). 



plot of the logarithmic derivative of the data (Fig. [6jb)) 
shows that the form of the curves close to the fixed point 
may correspond, in the large system size limit, to have 
the same exponent as the one of the surface scaling. 

Note that the (3 ~ value for the free surface case 
means that probably in the present model there are not 
significant difference between the split bottom and split 
top bottom geometry of double size. To confirm this fact, 
in Fig. [7| I show the results for the split bottom, and 
split top-bottom configuration with H = 80, where it 
can be seen than there are not any significant differences 




FIG. 5: (Color online) (a)Full profile of the width of the 
shear band across the sample for samples of different size. 
(b)Evolution of the value of w at the surface with system 
size. The error bar is smaller than symbol size. The contin- 
uous line is the best fitting, with an exponent 0.68. Dotted 
lines have slope 0.60 and 0.75, for comparison. 



between the two. This results was not obvious from the 
beginning, since the symmetry of the split top bottom ge- 
ometry with respect to the horizontal middle plane exists 
only on average, and not instantaneously. The present 
finding agrees with the observation in Ref. 6 of a profile 
that hits the surface with zero derivative (in that work 
however, it was not checked that the full profile was iden- 
tical to half profile of a split top bottom case with double 
size). In Ref.— it was also observed that the full profile 
of the w(x) could be nicely fitted with a quarter of a 
circle profile. An attempt to do this in the present case 
shows (see dashed line in Fig. [6]) apparent differences 
between this function and the numerical data that seem 
to be much larger than those reported in 6 -. I will come 
to this attempt of fitting in the next section. 

In addition to the plot of the width of the shear band, 
it is interesting to evaluate the full profile of the velocity 
at some typical vertical position, i.e, to consider the ve- 
locity v(y) at some fixed value of x. Previous numerical 
and experimental results show that to a good approxima- 
tion this profile can be fitted by an error function. The 
experimental data do not have sufficient precision so as 
to observe systematic differences with respect to a perfect 
error function. On the other hand, the numerical simu- 
lations in 6 show profiles that go above the error function 
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FIG. 7: (Color online) Comparison between the profile ob- 
tained for a split bottom configuration (free end case), and a 
split top-bottom configuration (fixed end case), for H = 80. 
No systematic differences can be observed within the numer- 
ical errors. 



FIG. 6: (Color online) (a)The results of the previous fig- 
ure scaled using an exponent a — 0.68. The fitting using 
Gaussian signals and a quarter of circle profile are shown. 
(b)Logarithmic derivative of the data in (a), and the corre- 
sponding curve from Gaussian signals and for a circle. 



for horizontal positions far away of the center of the shear 
band. It is not totally clear however if this difference will 
persist in infinitely long runs, as the difference slightly 
reduces as the simulation time increases. 



Given the velocity profile v(y) at some fixed x position, 
its y derivative will show a Gaussian form. This quantity, 
noted a(y) is presented in Fig. [5] for the plastic model, 
and for DPs. Data at the surface of the system, and 
in some particular positions in the bulk are included to 
show that there are no systematic differences between 
surface and bulk. The y axis is scaled in each case by 
the corresponding value of w(x), and the normalization 
cr(0) = 1 is used, in such a way that the best Gaussian 
fitting is given in all cases by the function exp(— x 2 /2). 
The results for DPs fall clearly below the Gaussian fitting 
for sufficiently large values of the argument, namely in 
the tail of the distributions. This is a well known result. 
Results for the full model indicate a very good Gaussian 
fitting for \y\ < 2w(x). Data seem to indicate a slight 
deviation in excess for larger values of y. However, these 
result require more massive simulations to be confirmed. 
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FIG. 8: (Color online) Derivative of the velocity profile a(y) = 
dv(y) /dy for the directed polymer (full symbols) and the plas- 
tic mesoscopic model (open symbols). Data at x = -ff/4, 
x = H/2, x = 3i//4, an x = H are put all together by scaling 
the y axis with the corresponding values of the width w. The 
deviation of DPs from a Gaussian (in defect for large values 
of x) is clearly visible. The quality of the data for the plastic 
model is not sufficient to make a statement about deviations 
with respect to a Gaussian. 



C. Plastic mesoscopic model: three dimensions 

The results of the previous section indicate that the 
two dimensional plastic mesoscopic model has similari- 
ties and differences with the fluctuating shear line, or the 
DP model, and with results of other kind of simulations 6 . 
In particular, the result obtained for the exponent a = 
0.68 ± 0.02 is remarkably consistent with the analytical 
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2/3 exponent of DPs. Experimentally, the obtained expo- 
nent is consistent with a 2/3 value, but to my knowledge 
an independent determination has not been attempted, 
and probably values within ±0.1 of 2/3 will be equally 
consistent with the experiments. 

However, an important point of discrepancy between 
the present results and atomistic numerical simulations 
in Ref£ is the form of a quarter of a circle profile found 
for w(x) in£. The results of the previous Section (see Fig. 

seem to depart from this behaviors in a larger extent 
than the results in£. 

One of the possible reasons for this discrepancy is the 
fact that the results m£ and experiments in general are 
fully three dimensional, while all the results presented 
up to now here are two dimensional. It is important to 
emphasize that although in the geometry studied there 
is translation symmetry on average along the z direction, 
the stochastic nature of the problem implies that fluctua- 
tions are not invariant along this direction. Actually, the 
setting up of a model that can be implemented in three 
dimensions was one of the motivations of the present re- 
search. 

A full implementation of the mesoscopic plasticity 
model in three dimensions requires the use of a three 
component displacement field Ui(x, y, z), i = x, y, z. This 
presents some technical challenges that I have not over- 
come yet. However, a partial implementation of the three 
dimensional case is possible, and the results obtained are 
already very interesting. 

The idea is to keep working with a single component 
displacement field u which is the displacement along z 
direction, keeping u x — u y = 0. The model can thus be 
described as a collection of planes in which a set of equiv- 
alent two dimensional problems are defined, and where 
the plastic thresholds are chosen independently for each 
plane. The fundamental variable is then some u{x,y)i 
where index i labels the plane along the z direction. 
Planes are connected to nearest neighbor ones by purely 
elastic forces, in such a way that the total energy of the 
system is now 



F = E ( F ' + y / dxd y^ 



H+l 



(10) 



where F l is the previous two dimensional free energy for 
plane i. Periodic boundary condition are applied along 
z. The size of the system along z should be chosen large 
enough. To keep a consistent implementation, this size is 
always chosen to be equal to the total thickness H of the 
system. The coupling constant k z along z is chosen to 
be equal to the linear elastic constant within the plane. 
Note that no plastic threshold is implemented for the 
springs along z. This version of the model was simulated 
in exactly the same manner as the two dimensional case. 
The results for the width function w(x) for systems of 
different sizes are shown in Figs. [S] and [11)1 Data for dif- 
ferent values of H can again be scaled by plotting w / H a 
vs w/H. The fitting of a gives a = 0.60 ± 0.02. 
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FIG. 9: (Color online) (a)Three dimensional data for different 
values of H. (b)Evolution of the value of w at the surface 
with system size. The error bar is smaller than symbol size. 
The continuous line is the best fitting, with an exponent 0.60. 
Dotted lines have slope 0.55 and 0.65, for comparison. 



Two important differences appear between the three 
dimensional and the previous two dimensional simula- 
tions. First of all, the scaling of data for systems of differ- 
ent thickness is obtained using an exponent smaller than 
that of the two dimensional case, namely a 3D ~ 0.60, 
a 2D ~ 0.68. The precision of the present simulations 
allows to claim that this difference is significant. In my 
view other numerical simulations as those in 4 or experi- 
ments have not reached yet the point to tell the growing 
exponent with a precision of ± 0.1. I think a value close 
to 0.6 is consistent with the precision of accessible exper- 
imental results. More experimental work is needed con- 
cerning this point. The second difference between two 
and three dimensions is the overall shape of the w(x) 
function. The form of this function in 3D is much closer 
to a quarter of circle profile (see Fig. [TO)) . This may be an 
interesting observation when comparing with the results 
in£, for instance. However, I have to mention that the 
fitting to this functional form does not seem to improve 
as the system size is increased, systematic differences ap- 
parently exist. 

Related to this behavior, there is an interesting point 
to emphasize. Given a set of curves w(x) for different 
thicknesses H, if the scaling hypothesis holds we have 



w{x) = w {y/H)H a 



(11) 
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FIG. 10: (Color online) (a)The data of the previous figure 
corresponding to the largest system size (H = 120) are shown 
together with a fitting using a Gaussian signals with exponent 
0.60 and a quarter of circle profile, (b) Logarithmic derivative 
of the data in (a), and the corresponding curves from Gaussian 
signals and for a circle. 

where wq is an TJ-independent function. If we assume 
wo behaves as a power law (with exponent a') for small 
values of its argument, we have 

w{x -> 0) - (x/H) a 'H a = x a ' H a - a ' (12) 

Note that unless a = a' , the behavior of w close to the 
origin will depend on the value of H, no matter how 
large this is. In particular, w(x — > 0) will increase with- 
out limit (if a > a') or decrease to zero (if a < a') as H 
increases. This would be in my view a very strange be- 
havior, although I have no rigorous argument against it. 
The DP curves shown in Figs. [2] and [3] are in fact com- 
patible with a 2/3 power close to the fixed point, and for 
the two dimensional plastic mesoscopic model a value of 
a' = a seems to be compatible with data in Figs. and 
[5J The situation is less clear in the three dimensional 
case. The value of a was found to be around 0.60. How- 
ever, the value of a' from Fig. [TU]seems to be consistently 
closer to 0.5 (which, by the way, is why the quarter of 
circle profile is accurate). It is not clear if the two ex- 
ponents will show a tendency to become closer in larger 
system size simulations of if they will remain different. 



IV. SUMMARY AND CONCLUSIONS 

In this paper I presented results obtained for the 'split- 
bottom' geometry using a mesoscopic model for the plas- 
ticity of an amorphous material. Two and three dimen- 
sional results were reported. In particular, the exponent 
for the width w of the shear band at the surface of the 
system for two dimensions was found to be 0.68 ± 0.02, 
comparable with the analytical 2/3 exponent obtained 
for directed polymers. In three dimensions the exponent 
was found to be somewhat slower, namely 0.60 ± 0.02. 
I want to recall however that the implementation of the 
three dimensional case was only partial, and a full mod- 
eling would be necessary to verify if this value remains 
valid in the full implementation. In two dimensions, the 
full profile of the width in bulk w(x) was found to be well 
described by a random Gaussian signal, and in particu- 
lar it was found that no systematic differences appear 
between a system with a free top surface and system 
with a split surface at the top. This was not an obvious 
point from the beginning as the reference system for this 
study, namely the directed polymer, shows systematic 
differences between the two cases that I have addressed 
and described in terms of the theory of random Gaussian 
signals. In three dimensions, it was observed that the 
form of w(x) can be fitted reasonably well by a quarter 
of circle profile, a result that has been obtained in pre- 
vious work in a related model^. A caveat about whether 
this scaling can continue to be valid for arbitrarily large 
systems was put forward. 

The number of experimental results that the present 
model describes, added to the very simple ingredients 
on which it is constructed, places it as a good candi- 
date in which the deep origin of the regularities emerg- 
ing from experiments can be elucidated. One promising 
route would be trying to generate continuous differen- 
tial equations that describe the evolution of probability 
densities of the stochastic model. This is an interesting 
prospect for future work. In addition, I think that the 
results presented can be a motivation for experimental- 
ists to obtain more accurate data of the quantities at 
play, to investigate how the values presented here fit the 
experiments. 
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